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ABSTRACT 


/—-^A  sharper  Sudakov-type  lower  bound  for  the  lower  confidence 
limit  on  the  reliability  of  a  series  system  than  the  one  given  in 
Harris  and  Soms  (1980)  is  obtained.  Numerical  examples,  coverage 
probabilities  and  the  listings  of  the  short  FORTRAN  programs  used 


are  also  provided. 
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SIGNIFICANCE  AND  EXPLANATION 


Serial  Systems  arise  naturally  in  practice  in  engineering  and 
physics.  Therefore  it  is  of  substantial  significance  to  be  able  to 
efficiently  utilize  data  obtained  on  individual  components  for  the 
purpose  of  obtaining  an  overall  assessment  of  the  reliability  of 
the  system.  The  methods  of  this  paper  may  be  employed  for  this 
purpose.  As  an  example,  suppose  there  are  five  subsystems,  the 
sample  sizes  are  (10,  15,  20,  25,  30)  and  the  observed  failures 
(1,  2,  1,  1,  2).  Then  the  methods  of  this  paper  give  that  the 
optimum  95%  lower  confidence  bound  on  the  system  reliability  lies 
between  .3871  and  .3934. 


- - " 

The  responsibility  for  the  wording  and  views  expressed  in  this 
descriptive  summary  lies  with  MRC,  and  not  with  the  authors  of 
this  report. 


IMPROVED  SUDAKOV-TYPE  BOUNDS  FOR 
OPTIMAL  CONFIDENCE  LIMITS  ON  THE  RELIABILITY 
OF  SERIES  SYSTEMS 

★  ★  * 
Bernard  Harris  and  Andrew  P.  Soms 


1.  The  Improved  Sudakov-Type  Bound 

We  will  adhere  to  the  notation  of  Harris  and  Soms  (1980)  and 

we  refer  the  reader  to  this  report  for  background  material  and 

additional  references.  Let  V.,  1  <  i  <  k,  be  independent  binomial 
random  variables  with  parameters  n^  and  p.  (the  success  probabil¬ 
ity)  and  let  X.=n^-V^,  *  1  -  ,  1  <  i  <  k with  n^<n2<...<nk 

here  and  throughout.  It  is  desired  to  obtain  a  "good"  lower 

k 

confidence  limit  on  II  p . ,  the  reliability  of  a  series  system 

1-1  1 

with  Independent  components.  Let  g(x)  ■  g(x.j  ,x2, . . .  .x^)  be  a 
monotonlcally  decreasing  or  Increasing  (In  each  component)  order¬ 
ing  function  for  real  (not  necessarily  Integral)  x.,  0<x^<n<t 
with  large  or  small  values,  respectively,  of  g(x)  being  best. 

Let  Tj  > r2 >  . . .  > rs  be  the  ordered  values  of  g(x)  In  the  decreas¬ 
ing  case  and  r^<r2<...<rs  In  the  1  ncreasl ng  and  1  et 
*1  "  {*1  *  ri}»  1  ■  1 ,2,...,s.-  Then  (Aj  ,A2,.. .  ,A$)  Is  a 

monotonic  partition,  l.e.,  (0,0, . . . ,0) e  A1 ,  ( n ^  ,n2, . . . ,n. ) e  A 
and  If  3t1  «  (x^  ,...^xlk),  x2  *  (x?1 , . . .  ,x2k)  with  *u<*21* 
1*1,2,.  ..,k,  then  Xj  eA|  Implies  x2  e  Aj ,  j  >  1 . 

— -  -  • 
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k  ,nil  v(  nrvf 
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and  for  1  <  n  < s-1 ,  let 


•„  -  <*»{.vi|  i  .  i  ft;.;?)  *  4  <2-j> 

M»1  1<n  Xj:x.  eA.  J  / 

J  J  * 

and  a  *  0.  Each  such  partition  may  be  identified  with  a  function 
defined  on  the  set  of  sample  outcomes  by  defining  the  ordering 
function  g(x),  where 


g ( x )  *  n  if  x  e  An  ,  l<n<s  ; 


(2.3) 


obviously  g(x)  Inherits  the  monotonicity  properties  of  the 
partition. 

For  defining  Buehler's  (1957)  method  of  optimal  lower  confl- 
k 

dence  intervals  on  n  p,  we  assume  that  g(x)  has  been  redefined, 

1-1  1 

If  necessary,  as  in  (2.3).  Then  we  have,  from  Harris  and  Soms 
(1980), 

Theorem  1 .  Let  X  be  distributed  by  (2.1).  Then  ag(x)  *s  a  (^”a) 

lower  confidence  bound  for  II  p4.  If  b_/3\  Is  also  a  (1-a)  lower 

k  1-1  1  9VM 

confidence  bound  for  H  p.  which  Is  monotonical ly  decreasing  in 

1-1  1 

g(x),  then  b^<a^,  1  <  1  <  s . 

Me  now  let  g(x)  denote  the  original  ordering  function,  since 
this  Is  necessary  for  the  applications  below.  In  order  to  obtain 
bounds  for  a^j  j  we  must  assume  that  xQ  Is  such  that  for  each 
t  >  1,2,...,k,  the  equation 


»^t  *  9^  X0 ) 


(2.4) 


has  a  unique  solution  yt,  y^<nt,  where  1f,  r  -  1,2,...,t-l,  are 
Integers,  0<1r<yr.  Define  y*  by  g(0,0, , . .  ,0,y*  ,0 . 0)  •  g(xQ) 


K  X 


where  y^  Is  In  the  1th  position,  the  rest  of  the  arguments  being 
0.  Note  that  y,  =  y^.  A  sufficient  condition  for  (2.4)  to  hold 


is  that 


g(x)  is  strictly  monotonic  for  0  <  ,  j*l,2,...,k' 


Urn  g ( x )  <  g ( x-  )  for  g(x)  decreasing 

v\ 


(2.5) 


11m  g(x)  > g(x  )  for  g(x)  increasing, 

x+n.,  r=l  ,2 , . . . ,  k  . 

r  r 


Assume  now  that  (2.5)  is  satisfied  and  that  in  addition 


*r-’r  . 

"r-'r  -  nr+l  ' 


r«l  ,2 . k-1  , 


(2.6) 


where  1f  integral,  0<lriyr*  Let 

and  for  0  <  y  <  n,  real ,  define  u(n,y,a)  by  a  =  Iy^n  aj(n-y,y+l ) 
Then  it  was  shown  in  Harris  and  Soms  (1980)  that 


<  ,9(-o)  <  .«) 


(2.7) 


and  thus  If  y^  Is  an  Integer,  j  ■  u(n^,y^,a). 

Thus  a  conservative  procedure  Is  to  use  u(n^,y^,a)  for 
*9(x  )*  **  Is  generally  believed  that  this  Is  quite  conservative 

(see,  e.g.,  Mann,  Schafer  and  Slngpurwal la ,  1974).  It  Is  thus  of 
practical  Interest  to  inquire  whether  the  lower  bound  In  (2.7)  can 
be  tightened.  An  examination  of  the  proof  of  (2.7)  In  Harris  and 
Soms  (1980)  shows  this  to  be  the  case,  with  the  new  proof  being 
coincident  with  the  previous  one,  except  for  the  omission  of  the 
final  step.  We  have 


Theorem  2.  Under  the  same  assumptions  as  for  (2.7), 


u'fwyj.yj.o)  <  a,,^)  <  ^J"k  »(■•,. ty,].a) 


where  u '  (n^  .rig.y-j  Is  the  solution  In  a  of  the  equation 

^yl^  f  ni  1  ni"ii 

kS“P  I  .  P,  «1  1  k  (»2-y,.y2*U  -  «  • 

#1"*  V°  S-V  n  p, 


It  follows  from  the  proof  In  Harris  and  Soms  (1980)  that 


ufnj.yj.a)  <  u* (n^ ,n2,y1 ,y2,a)  . 

Thus  the  only  question  Is  how  great  the  Improvement  will  be 
In  using  u'.  The  numerical  examples  In  2.  show  that  this  can  be 
substantial. 

Remarks .  u'  can  be  calculated  quickly  and  efficiently  by  using  a 
short  FORTRAN  program.  The  listing  Is  given  In  the  Appendix  along 
with  the  listing  of  the  program  used  to  calculate  coverage  proba¬ 
bilities.  Two  ordering  functions  that  satisfy  (2.5)  and  (2.6)  are 
k  k 

g(x)  »  n  ((n.-x.)/n.)  If  g(x  )>0and  g(x)  *  7  x./h .  for  sufficiently  large 
1»1  1  1  1  0  1«1 
n1*  1  < 1 <  k  (see  Harris  and  Soms,  1980,  for  details). 


2.  Coverage  Probabll itfes  and  Numerical  Examples 


The  ordering  function  used  here  throughout  Is 
k 

g(x)  *  n  ( (n.-x. )/n. ) .  Table  1  gives  the  coverage  probabilities 
1-1_  111 

for  k«3,  n*  (5,7,10),  a  •  .10  and  selected  p  ■  (Pj^^K  for  both 
the  lower  LB  and  upper  UB  bounds  of  (2.7).  While  the  optimality 
property  of  Theorem  1  Implies  that  there  are  p  for  which  the 


coverage  probability  Is  less  than  .9,  It  does  not  seem  that,  there 


Coverages 


-P  -- 

UB 

LB 

(.95, .95,. 95) 

.9999 

.9999 

(.95, .95, .90) 

1.0000 

1.0000 

(.95, .95, .85) 

1.0000 

1.0000 

(.95, .95, .80) 

.9999 

.9999 

(.95, .95, .75) 

.9999 

.9999 

(.95, .95, .70) 

1.0000 

1 .0000 

(.95, .95, .65) 

.9535 

.9927 

(.95, .95, .60) 

.9737 

.9967 

(.95, .95, .55) 

.9869 

.9874 

(.95, .95, .50) 

.9940 

.9940 

Table  2  gives  some 

comparisons  of  LB, 

UB,  the  improved  lower 

bound  LBI  and  the  true 

value  TV  for  k  =  2, 

a  =  .05  and  selected 

n  =  (n^.n^)  and  xQ  =  (x 

o i »x02 ) •  LBI  9ives 

substantial  improvement 

when  n^  is  small  compared  to  n^  and  little 

or  none  if  n^  is 

approximately  the  same 

as  n2-  The  TV  for 

n  =  (5,5)  and  x  =  (1.1) 

agrees  with  the  value  in  Lipow  and  Riley  (1959). 


2.  Comparison  of  Bounds  and  True  Value  for  k  =  2  and  ot  «  .05 


n 

X  1 

o 

LB 

LBI 

TV 

UB 

(5,5) 

(1.1) 

.2166 

.2166 

.2776 

.3426 

(5,10) 

(1,1) 

.2761 

.3317 

.3317 

.3426 

(5,10) 

(2.3) 

.0859 

.1521 

.1529 

.1893 

(10,20) 

(1,2) 

.5037 

.5675 

.5691 

.6058 

(10,20) 

(4.5) 

.1851 

.2137 

.2137 

.2224 
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Table  3  gives  LB,  UB  and  LBI  for  k  *  5,  a  s  .05  and  selected 
n  and  x. 


3.  Comparison 

of  Bounds  for 

k  *  5  and  a  * 

.05 

n 

*o 

LB 

LBI 

UB 

(10.10,10,10,10) 

(1,2, 0,1,1) 

.2893 

.2893 

.3035 

(10,15,20,25,30) 

(1.2, 1,1, 2) 

.3599 

.3871 

.3934 

(10,15,20,25,30) 

(1.1. 2, 3, 4) 

.2838 

.3017 

.3035 

(10,15,20,25,30) 

(2, 3, 4, 4, 4) 

.1319 

.1478 

.1500 
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Appendi x 


The  listings  of  the  coverage  and  improved  boun  j  FORTRAN 
programs  are  given  below. 


N’S  SAMPLE  SIZES  FROM  SMALLEST  TO  LARGEST,  P'S  SUCCESS  PROBABILITIES 
K  NUMBER  OF  SUBSYSTEMS,  1-ALPHA  CONFIDENCE  LEVEL 

DIMENSION  P ( 1 0 ) ,N(20)  ,NY(10)  ,  PL IMU ( 1 0 ) , P I ( 1 0 , 41 ) , LL ( 1 0 ) , 1 1 ( 1 0 ) 

100  FORMAT  (  ) 

EQUIVALENCE  ( 1 1 , 1 1 ( 1 ) ) , ( 1 2 , 1 1 ( 2 ) ) , ( 13 , 1 1 ( 3 ) ) , ( 14 , 1 1 ( 4 ) )  , 

1(15, 11(5)), (16, 11(6)), (17, 11(7)), (18, 11(8)), (19, 11(9)), (110, 11(10) 
2) 

1  READ  100, K, ALPHA 
IF  (K.EQ.O)  GO  TO  21 
DO  221  1=1 ,10 
LL ( I )  =  0 
221  N( I  )  =  0 

READ  100,(N(I),P(I),I=1,K) 

DO  343  111  =  1  , K 
DO  343  112=0, N(II1) 

CALL  MDBIN(II2,N( III ) ,P( III ) ,PT,PI( III ,112+1 ) ,IER) 

343  IF  (PT.LE. .0001  )  LL ( 1 1 1 )  =  LL ( 1 1 1 ) +  1 
PARAM= 1 . 

I PR0D2=1 

DO  11  1  =  1  ,  K 

IPR0D2=IPR0D2*N(I) 

11  PARAM=PARAM*P( I ) 

PRL=0 . 

PRU=0 . 

DO  122  I1=LL(1),N(1) 

DO  122  I 2=LL ( 2 ) ,N(2) 

DO  122  I3=LL(3) ,N(3) 

DO  122  I4=LL (4 ) , N( 4) 

DO  122  I 5=LL ( 5 ) ,N(5) 

DO  122  I6  =  LL(6)  ,N(6) 
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DO  122  1 7  =  LL ( 7  ) ,  N ( 7  ) 

DO  122  18  =  LL ( 8 ) ,N(8) 

DO  122  I 9=LL ( 9 ) ,N(9) 

DO  1 22  1 1 0  =  LL ( 1 0 ) , N ( 1 0  ) 

PROBT  =1 . 

DO  51 1  I T= 1  ,  K 
INT= 1 1 ( IT) 

511  PR08T=PR0BT*PI ( IT, INT+1  ) 

I PROD1 = 1 
DO  4  I J  =  1  ,  K 
I  N= 1 1 ( I J  ) 

4  I  PROD  1  =  I  PROD  1  *  I  N 

IF  ( IPROD1 . EQ . 0 )  GO  TO  15 

ITEMl=IPROD2-IPRODl 

Y 1  =  N ( 1 ) *  I TEM1 / FLOAT ( IPROD2) 

DO  5  I T  =  1  ,  K 

NY(IT)=N(IT)*ITEMl/IPROD2 
A=FLOAT(N( IT)-NY(IT)  ) 

R  =  FLOAT ( NY ( I T )  + 1  ) 

5  CALL  MDBETA(PARAM,A,B,PLIMU( IT)  ,IER) 

CALL  MDBETA( PARAM, N( 1 ) - Y 1 , Y1 +1 . , PL  I  ML , I ER ) 
IF  (PLIML.GE. ALPHA)  PRL  =  PRL  +  PROBT 
I FLAG=0 
DO  6  I T=  1  ,  K 

6  IF  (PLIMU(IT) .GE. ALPHA)  I FLAG= 1 
IF  ( 1  FLAG . EQ . 1 . )  PRU  =  PRU  +  PROBT 
GO  TO  122 

15  PRL=PRL+PROBT 
PRU=PRU+PROBT 
122  CONTINUE 

PRINT  100, K, ALPHA, (N( I) ,P( I) ,1=1 ,K) 

PRINT  1 00 , PRU , PRL 
GO  TO  1 
21  STOP 
END 
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C  SERIES  WITH  UNLIKE  COMPONENTS 

C  MDBIN,  MDBETA ,  MDBETI  IMSL  BINOMIAL,  BETA,  INVERSE  BETA  ROUTINES, 

C  N'S  SAMPLE  SIZES  FROM  SMALLEST  TO  LARGEST,  X ' S  FAILURES,  NINT  NUMBER 
C  OF  P  VALUES  CONSIDERED  FOR  SUP,  EPS  MAXIMUM  ERROR  IN  SOLUTION  FOR 
C  REFINED  BOUND,  NO  X'S  EQUAL  TO  N’S 
REAL  ITEM1 ,ITEM2,ITEM3,ITEMT 
DIMENSION  N ( 50 ) , X ( 50 ) 

100  FORMAT  (  ) 

DATA  EPS/. 001/ 

1  READ  100,  K, ALPHA, NINT 
IF  (K.EQ.O)  GO  TO  99 

READ  100,  (N(I) ,X(I),I  =  1 , K ) 

52  PRINT  100, 'DATA' , (N( I ) ,X( I ) ,1  =  1  ,K) 

PRINT  100,  '  ALPHA= ' .ALPHA 
I TEM1  =  1  . 

I TEM2  =  1  . 

DO  2  1  =  1  ,K 
ITEM1  =  ITEM1 *N( I  ) 

2  I TEM2= I TEM2* ( N ( I  )-X(  I ) ) 

ITEM3=N( 1 )*( ITEM1 -ITEM2) 

Y 1 = I TEM3/ I TEM1 

NY  1  =  Y  1 

CALL  MDBETI ( ALPHA, N( 1 )-Yl  ,Y1+1 . ,A,IER) 

C  COMPUTATION  OF  UPPER  BOUND 
BU  =  1  . 

DO  3  1  =  1  ,K 

ITEMT=N( I )*( ITEM1 -ITEM2) 

NYT=ITEMT/ ITEM1 

CALL  MDBETI (ALPHA, FLOAT (N( I )-NYT)  ,NYT  +  1 .  ,BUT,IER) 

3  IF  (BUT.LT.BU)  BU=BUT 

PRINT  100, 'L-M  LOWER  BOUND  =',A 
PRINT  100,  'UPPER  BOUND  =' ,BU 
C  COMPUTATION  OF  IMPROVED  LOWER  BOUND 
TEML=A 
TEMU=RU 

9  TEM= ( TEML+TEMU ) /2 . 

DELTA= ( 1 -TEM)/NI NT 
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SUP=0. 

DO  42  1=1 ,NINT 

PI =TEM+( I -1 ) *DELTA 

P2=TEM/P1 

Q1 =1 -PI 

PR0B=0 . 

DO  33  K1 =0 , N Y 1 

Y2=N(2)-( ITEM2/ITEM1 )*( N( 1 )*N( 2 )/ ( N( 1 ) -K1 ) 

CALL  MDBIN(K1 ,N(1 ) ,Q1 ,PS1 ,PK1 ,IER) 

CALL  MDBETA(P2,N(2)-Y2,Y2  +  1 .  ,  PP2 , 1 ER) 

33  PR0B=PR0B+PK1*PP2 

IF  (PROB.GT.SUP)  SUP=PROB 
42  CONTINUE 

IF  (SUP. LE. ALPHA)  TEML=TEM 
IF  (SUP. GT. ALPHA)  TEMU=TEM 
IF  (ABS(TEMU-TEML) .LE.EPS)  GO  TO  10 

GO  TO  9  | 

10  A=TEML 

PRINT  100, ' IMPROVED  LOWER  BOUND  =  '  ,A 
GO  TO  1 
99  STOP 
END 
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